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In the ferromagnetic insulator with Dzyaloshinskii-Moriya interaction, we theoretically predict and numeri- 
cally verify a topological magnon insulator, where the charge-free magnon is topologically protected to transport 
along the edge while it is insulating in the bulk. Within the bulk band gaps, edge states form a connected loop as 
a 4.7r- or 8-7r-period Mobius strip in the wave vector space. As a consequence, the chiral energy current traveling 
along the corresponding edge is topologically protected from defects or disorders. Using the nonequilibrium 
Green's function method, we demonstrate that the energy current carried by magnons with energy in the bulk 
gap localizes at edges and prefers to travel along one edge in only one direction at nonequilibrium steady state. 
Our prediction about topological magnon insulator could be observed at a wide energy range in the thin film of 
the insulating ferromagnet, such as LU2V2O7. 

PACS numbers: 85.75.-d, 75.30.Ds, 75.47.-m, 75.70.Ak 



Topological insulator, as a novel state of quantum matter, is 
characterized by an insulating bulk band gap and conducting 
gapless edge/surface states protected by symmetries [1, 2]. It 
has been theoretically predicted and experimentally observed 
in a variety of systems and becomes a hot spot because of its 
theoretical importance in condensed matter physics and wide 
potential applications in dissipationless spin-based electron- 
ics (spintronics) [3]. However, due to the fact that the spin 
transport in topological insulators is carried by electrons, dis- 
sipations can not be really avoided. 

Magnon Hall effect, as a consequence of the 
Dzyaloshinskii-Moriya (DM) interaction [4, 5] which 
plays a role of vector potential similar to the Lorentz force, 
has been predicted and observed in magnetic insulators 
[6-8]. Compared with spin current, where the dissipation 
from Joule heating is still inevitable due to the electronic 
carriers, the magnon Hall effect is more promising in device 
applications because of the long-range coherence of charge- 
free spin wave [9-11]. Magnons are collective excitations 
of localized spins in a crystal lattice and can be viewed as 
quantized quasiparticles of spin waves. Recently, magnon 
excitation [12, 13], localization [14] and interference [15] 
have been experimentally realized thanks to the achievements 
in deposition and nanopatterning techniques. The technical 
advancements offer the perspective of various magnonic 
devices, and a new discipline - magnonics - has emerged 
and growing exponentially [16-19]. The charge-free property 
of magnon makes it promising to achieve dissipationless 
transport and control in insulating magnets. 

Therefore, it will be of great general interest for both the- 
orists and experimentalists that we find in this work a new 
intriguing quantum state that magnon while insulated in the 
bulk, can nondissipatively transport along edges/surfaces in 
the absence of backscattering from defects and disorders due 
to the nontrivial topology of magnon's band structures. We 
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FIG. 1: (color online), (a) Pyrochlore crystallographic structure of 
the sublattice of magnetic atoms V of LU2V2O7; (b) Two tetrahe- 
drons in the pyrochlore lattice, where DM vectors on bonds 1-3, 1-2, 
and 2-4 are indicated by orange arrows; (c) Schematic magnetic flux 
due to DM interaction in the [111] plane of the pyrochlore lattice, 
i.e., a kagome lattice. The coupling of two sites along the arrows is 
(J + iD)S, and the opposite direction corresponds to (J — iD)S; (d) 
The quasi-one-dimensional strip kagome lattice. The area enclosed 
by the dotted line can be regarded as a center which is connected to 
two semi-infinite leads in equilibrium at temperatures Tl (left) and 
Tr (right). The two big arrows along edges schematically depict the 
energy flow magnitudes and directions in the lattice when Tl > Tr. 
The width of the strip is W = 5, and the number of atoms in each 
unit cell is 6W — 1. 



name this novel state topological magnon insulator (TMI) and 
believe that due to the robust dissipationless magnon trans- 
port, the TMI in insulating magnets could provide widely po- 
tential applications in nondissipative magnonics and micro- 
spintronics. 

The magnon Hall effect was experimentally observed in in- 
sulating ferromagnet LU2V2O7 [7] with a pyrochlore lattice, 
in which the magnetic atom vanadium has a corner-sharing- 
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tetrahedra sublattice, that is, a stacking of alternating kagome 
and triangular lattices along the [111] direction, as shown in 
Fig. 1 (a). To study magnon transport in the ferromagnetic 
insulator, the Hamiltonian can be written as [4, 5, 20]: 
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where S n is the spin angular momentum at site n; — J denotes 
the nearest-neighbor coupling; D mn is the DM interaction be- 
tween site m and n; the last term comes from the Zeeman 
effect under an external field H. 

As shown in Fig. 1 (b), in a single tetrahedron, the DM 
vector is perpendicular to the corresponding bond and paral- 
lel to the surface of the surrounding cube [7, 21, 22]. Since 
the component of D mn perpendicular to z = H/H does not 
contribute to the Hamiltonian (1) up to quadratic order in 5S 
[7], we only retain D mn = D mn ■ z. When applying a mag- 
netic field along z = [111] direction, all the projections of 
the DM interaction between inter-layer sites m and n are zero 
(Di3 = D23 = -D43 = 0); and all the ones for inner-layer 
sites are nonzero and equal (D12 = -D24 = -D41 7^ 0). There- 
fore, with the magnetic field along [111] direction, magnon 
Hall effect only comes from the noncancellation of different 
types of loops of DM interaction in the unit cell of the kagome 
lattice, as shown in Fig. 1 (c). Thus the kagome sublattice 
structure plays a key role for the presence of magnon Hall ef- 
fect in L.U2V2O7. In the following we first discuss a general 
kagome lattice with DM interaction; latter we will incorporate 
actual parameters of a thin film of L.U2V2O7 with a kagome 
layer vanadium sublattice structure. 

Using the spin-wave (Holstein-Primakoff) transformation 
[23, 24], one can straightforwardly obtain the quadratic spin- 
wave Hamiltonian as [25] : 



H = 



+ E 



(2) 



where bVb denotes the operator that increases/decreases the 
spin component along z = [111] which is also the direction 
of the applied magnetic field, H mn = H* m = ( J ± iD)S, 
m n, H mm = JSAI n + hS, M n is the number of 
nearest neighbors of site n. The ground-state energy reads 
E = -JS 2 J2„ M «/2 - NhS with N the total number of 
sites and h = gfigH. We set D mn = D and the kagome 
lattice with DM interaction is shown in Fig. 1 (c), where the 
coupling between two sites along the direction of that arrow 
corresponds to (J + iD)S, and the opposite direction cor- 
responds to (J — iD)S [25]. Due to the different types of 
loops in a unit cell of the kagome lattice, the DM interaction 
avoids cancellation thus induces the Hall effect [6]. As a con- 
sequence, the DM interaction acts as a vector potential for the 
propagation of magnons similar to the magnetic field for the 
propagation of electrons. In the following calculations, we 
use dimensionless units and set S = 1/2, and J = 1 without 
loss of generality. The magnetic field decides the direction of 
spins at the ground state, and the induced Zeeman effect term 
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FIG. 2: (color online). The upper part of the dispersion relation of 
kagome lattice strip in the range of k x = ~ 4n/a (a is the lattice 
constant), with (a) D/J = 0.1, and (b) zero DM interaction. 



shifts the dispersion relation. We set magnetic field H = + 
in the first part of theoretical calculations, and will input finite 
H in the later part for real-material calculations. 

With the bosonized Hamiltonian (2), we are able to cal- 
culate the magnon's band structure as well as the Berry cur- 
vature and Chern number for each band [26]. To investigate 
the edge magnon transport, we carry out the calculation and 
analysis based on the strips of two-dimensional kagome lat- 
tice structure as shown in Fig. 1 (d). The dispersion relation 
of magnons (e vs k x ) in the infinitely long strip with DM inter- 
action is displayed in Fig. 2 (a). Since there are three sites in 
each kagome unit cell, there are three energy bands separated 
by two energy gaps. With increasing strip width, more modes 
appear in the energy bands; the edge states will be gradually 
fixed and independent of size [27]. When the system width is 
finite, the states in two edges could couple together to open an 
energy gap, which exponentially decays with width increas- 
ing [25, 28]. When the strip width increases to infinity, the 
edges are separated too far to interact with each other; thus 
they could have degeneracy in dispersions, as shown in Fig. 2. 
In our calculation, the edge states are almost invariable when 
the width W > 20 [25] and we use width W = 80 in all the 
latter numerical calculations. Because of the DM interaction, 
two edge states within both energy gaps are twisted so that 
for each state eiixja — k x ) 7^ e(tt /a + k x ), and they cross 
at k x = it /a in the first Brillouin zone, as shown in Fig. 2 
(a). As a consequence, the two edge modes form a continu- 
ous state with a period of 4n, which can not be disturbed to 
open a gap by weak disorders so that edge modes are topolog- 
ically protected. However, when DM interaction is zero, the 
two edge states are easy to be perturbed to separate and open 
a gap, because they do not cross with each other although they 
degenerate as shown in Fig. 2 (b). 

In the lower bulk band gap, we find that four edge states 
will contribute to magnon transport within the energy gap, as 
shown in Fig. 3 (a). If D/J « 0.4 or larger, in both energy 
gaps there are four edge states. In a period of 871", any two 
of the four edge states have degeneracies and cross each other 
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FIG. 3: (color online). The dispersion relation of edge modes with 
non-zero DM interactions, (a) and (b) are dispersion relations in 
range of k x = ~ 8n/a for D/J = 0.1 and 0.4, respectively, 
(c) A conventional cylinder strip with two boundaries, of which each 
has a period of 2ir; (d) A Mobius strip which only has one bound- 
ary of 4-7T period ; (e) A looped Mobius strip which only has one 
boundary of 8ir period. 



at different points in the momentum space. All the four edge 
states form a continuous state with a period of 8ir to transport 
magnons along two edges of the strip. We can understand the 
topology of the edge states as follows. With zero DM inter- 
action, the edge states are similar to the two boundaries of 
the conventional cylindric strip as shown in Fig. 3 (c), both of 
which have a period of 2ir in Brillouin zone and transport sep- 
arately along two edges. Due to the nonzero DM interaction, 
two edge states are twisted so that they cross each other and 
go into the other energy band after 2ir in momentum space, 
thus form a closed loop with a period of 47r. These edge states 
are similar to the one-sided Mobius strip with only one bound- 
ary, as shown in Fig. 3 (d), where a line drawn starting from a 
point at the boundary will meet back at the "other side" after 
a circle of 2tt, then go back to the original point after a whole 
period of 47r. The two edge states form one closed loop wind- 
ing the bulk energy gap between the two bands, which are 
thus topologically protected from distortions. At larger DM 
interaction, four edge states contribute to the transport in the 
bulk gap, cross each other, and connect to form a closed 8ir- 
period loop which can be interpreted as the only one boundary 
of a looped Mobius strip as shown in Fig. 3 (e). This looped 
Mobius strip also has only one boundary winding around the 
strip surface, thereby having the same topology as that of the 
conventional Mobius strip shown in Fig. 3 (d). 

The edge state transport is related to the band topology 
characterized by Chern numbers of the bulk states [27, 29- 
31]. When the DM interaction is absent, all the Chern num- 
bers of three bands are zero so that there is no magnon Hall ef- 
fect [25]. Accordingly, the winding numbers of edge states are 
zero as well, thus they are not topologically protected. With 
nonzero DM interaction, the Chern numbers of the lowest and 
highest energy bands become ±1; the one of the middle en- 
ergy band is still zero [25]. According to the relation between 
the Chern number and the winding number [27], the winding 
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FIG. 4: (color online). The local energy current and density of state 
for edge magnon transport at equilibrium, (a) The uniform kagome 
lattice, (b) The lattice with a defect at the upmost site of the left 
fifth column. The red arrows, the blue dots, and the small black dots 
correspond to the local energy current, the local density of magnon, 
and atom sites, respectively. The color of the arrow indicates the 
magnitude of the local current. Parameters are e — 1.5, Tl = Tr = 
1.0, D/J = 0.1, lattice spacing a = 1. 



numbers of edge states in both bulk gaps have the same value 
of 1 or —1, which is consistent with the only one closed loop 
winding the bulk gap regardless of the period of 47r, 8tt or 
others. 

To intuitively illustrate the magnon transport by edge states, 
we choose some unit cells of the kagome lattice strip as a cen- 
ter region and set the rest two semi-infinite lattices in equilib- 
rium at temperatures Tj, and Tr to study the energy transport 
by magnons as shown in Fig. 1 (d). Fig. 4 shows the edge 
state magnon transport in the bulk gap at a fixed magnon en- 
ergy e = 1.5 in thermal equilibrium. The forward (left-to- 
right) thermal current carried by magnons travels along one 
edge, and the backward (right-to-left) current with the same 
magnitude transports along the other one, as shown in Fig. 4 
(a). Near both edges the local magnon currents form many 
chiral vortices due to the nonzero DM interaction. Moreover, 
both the current and the magnon density of states are symmet- 
rically localized at two edges. They exponentially decay from 
the edge to the center with oscillations [25], which is simi- 
lar to the properties of the localized edge phonon modes [32] 
and electron transport in graphene [33]. If the DM interaction 
reverses (D — > —D), all the currents reverse to the opposite 
direction. 

When a defect is present at one edge, the current take a 
detour around it and transports ahead without backscattering, 
as illustrated in Fig. 4 (b). Although the defect dramatically 
affects the local density of magnons and destroys the local 
current vortex, the global currents along two edges keep intact 
compared to those of the uniform lattice, and their summation 
keeps zero since the net transport vanishes at equilibrium. 

For the nonequilibrium case that two leads are held at dif- 
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FIG. 5: (color online). The local energy current and density of state 
for edge magnon transport at nonequilibrium. The uniform kagome 
lattice are with (a) T L = 1.2, Tr = 0.8, and (b) T L = 0.8, Tr = 
1.2. The lattice with a defect at the upmost site of the left fifth column 
are with (c) T L = 1.2, Tr = 0.8, and (d) T L = 0.8, Tr = 1.2. 
Other parameters are e = 1.5, D/ J — 0.1, a = 1, W = 80. 



ferent temperatures, the magnon will prefer to flow along one 
edge in the steady state, see Fig. 5 (a) and (b). For positive 
DM interaction and Tj, > Tr, the forward thermal current 
transported along the lower edge is larger than the backward 
one along the upper edge. The summation of forward and 
backward currents along two corresponding edges satisfies the 
Landauer formula. If we swap the temperatures of the left and 
right end, the thermal current will prefer the other edge. If 
the defect is present, both the local density of magnon and 
the local thermal current will be affected. However, the one- 
way edge magnon transport can detour the defect and keep the 
same magnitude to flow ahead, see Fig. 5 (c) and (d). In this 
way, the global thermal current keeps invariant, which means 
that the edge state in the bulk gap is indeed topologically pro- 
tected from the lattice defect or disorder. 

With the help of nonzero DM interactions, the magnon chi- 
rally travels along the corresponding edge nondissipatively in 
the absence of backscattering, which is preserved from weak 
defects or disorders. We use the actual parameters of thin film 
of LU2V2O7, to investigate the energy range of magnons for 
observing the TMI. As shown in Fig. 6 (a), the energy cur- 



FIG. 6: (color online). The current density vs energy of magnon 
for uniform and edge-defect kagome lattices with the parameters of 
LU2V2O7. (a) The solid and dotted lines correspond to the thermal 
current in the bulk band gaps for uniform and edge-defect lattices, 
respectively, (b) The dispersion relation of the kagome lattice with 
the parameters of Lu 2 V 2 7 . D/J = -0.26, H = 1 T, T L = 21 K, 
and T R = 19 K. 



rent of magnon is not affected by defect or disorder in the 
energy range of [4.45, 5.98] U [8.79, 10.31] meV. These en- 
ergy ranges coincide with the bulk gaps in the energy spec- 
trum where the topological magnon edge states can be found 
[see Fig. 6 (b)]. Although certain distortions of edge states 
will occur as results of the defect or disorder, the energy cur- 
rent carried by magnons has no difference in the whole bulk 
energy gap. This indicates that the defect or disorder does 
not open a gap in the energy spectrum so that the topology 
of the magnon edge state is robust. According to the en- 
ergy ranges, the topological magnons have frequencies rang- 
ing from [1.08, 1.45] U [2.13, 2.49] THz. Applying different 
external magnetic fields will not change the main properties 
of magnons, but shift the corresponding dispersion relations 
[25], so that the frequency of topological magnons can be 
tuned flexibly with a wide range. Therefore, we expect that for 
the thin film of LU2V2O7 in a wide energy range of magnons, 
one can observe the TMI, which could also be highly likely 
for other magnetic insulators containing the kagome lattice 
structure. 

To realize magnonic devices as well as the predicted TMI, 
the excitation and detection of magnons is the major chal- 
lenge. Recent years have witnessed a fast development in ex- 
perimental techniques such as ferromagnetic resonance [15], 
pulse-inductive microwave magnetometer [34], time-resolved 
scanning Kerr microscopy [35], optical pump-probe tech- 
niques [36], as well as Brillouin light scattering (BLS) [37] 
which takes a special role since it allows the direct measure of 
dispersions and band structures. We could observe the topo- 
logical edge modes by using these techniques to measure the 
magnon dispersion relation and could also verify the TMI by 
detecting the magnon transport in the bulk band gap where 
magnons could be excited by non-thermal optical pulses [38] 
or spin-polarized current [39], so that we can avoid the ther- 
mal transport due to the bulk states thus extract only the edge 
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channels. We hope our theoretical predictions about TMI 
could open a new window into the application of nondissi- 
pative magnon transport, especially for the novel magnonic 
device design, which could also shed light on the information 
technology based on magnonics, and micro-spintronics. 
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